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Abstract 

We study equations with infinitely many derivatives. Equations 
of this type form a new class of equations in mathematical physics. 
These equations originally appeared in p-adic and later in fermionic 
string theories and their investigation is of much interest in mathe- 
matical physics and applications, in particular in cosmology. Differ- 
ential equation with infinite number of derivatives could be written as 
nonlinear integral equations. We perform numerical investigation of 
solutions of the equations. It is established that these equations have 
two different regimes of the solutions: interpolating and periodic. The 
critical value of the parameter q separating these regimes is found to 
be q 2 cr « 1.37. Convergence of iterative procedure for these equations 
is proved. 
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1 Introduction 



Normally, in mathematical physics one considers equation with finite number 
of derivatives. For such equations there are well developed methods of solu- 
tion for various boundary problems, see for example pQ. Recently in works on 
p-adic and then in real string theories a certain class of nonlinear equations 
which involve infinite number of derivatives is started to be explored j2]-[Zj- 
Such equations form an important new class of equations in mathematical 
physics. New methods to study uniqueness and existence of solution should 
be developed. It is not clear apriori how to pose boundary or initial value 
problems for these equations. 

An example of new equations has the form 

e aA V = (1) 

where A is the Laplace (or D'Alamber) operator, a is a real parameter, and k 
is nonnegative integer. This equation was originally studied in p-adic string 
theory. Soliton solutions to (JTJ) where considered in (21 Ej- It is interesting to 
study equation ((TJ) in the simplest case when <p = (p(t) depends only on one 
real variable t [Old- I n this case we can rewrite (JTJ in the form of integral 
equation 

(K V ){t) = V {t)\ (2) 

where 

e^V(t) = (^)(t) = -== / e-^<p(t')dt' (3) 

In recent paper [B] a kink solution found in j2] was confirmed and also oscil- 
latory solutions where found. 

In this paper the equations of the form (0) and more general equations are 
investigated. These equations describe dynamics of the scalar field with the 
lowest mass square (tachyon field) in fermionic string model jH] • We perform 
numerical investigation of solutions of these equations. It is established that 
these equations have two different regimes of the solutions: interpolating 
and periodic. The critical value of the parameter q separating these regimes 
is found to be q\ T ~ 1.37. Convergence of iterative procedure for these 
equations is proved. To construct numerical algorithm we essentially used 
object-oriented design, please see section l3~B1 for further details. 

This paper is organized as follows. In the next section we study the equa- 
tion taken from p-adic string model, we prove the convergence of iterative 
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procedure for this equation. In the next two sections we study equations 
describing the scalar field (tachyon field) in fermionic string. Equations for 
fermionic string generalize the equations for p-adic string. Finally, in the last 
section we provide physical details including actions for which corresponding 
equations of motion form the subject of this paper. 



2 Scalar Field Dynamics in P-adic String Model 

In this section we study the following nonlinear integral equation 

(Kjp)(t) = <p{ty, ( 4 ) 
where integral operator /C is defined by 

1 f°° 

(JC<p)(t) = -= / e-^r (p(f)d£ (5) 



Originally p in the rhs of (@J) is a prime number, although here it is not im- 
portant. This equation has several physical applications basically it describes 
dynamics of homogenous scalar field configurations in the p-adic string model 
El HJ We describe the physics behind in more detail in the section El The 
equation (J3J) for p = 3 was studied in j2J. It was numerically shown that 
it has solution which goes to =Fl on infinities. In this section we prove the 
convergence of iterative procedure used in [2] to construct numerical solution 
of (0J). In jH] it was shown that (@J) does not have monotonic solutions for 
even p. In the text below we discuss only the case p = 3 since it is the most 
illustrative for the p-adic string model and provides an approximation for 
fermionic string model. 

We consider equation (jU) for the case p = 3 

(/C^)(t)=^(t) 3 (6) 

We are searching for the solution which has constant asymptotic behavior on 
infinities. For the case of constant field equation (JHJ) takes the form 

fo = fl (7) 

It has three solutions 

= ^ 2) = 0, ^ 3) = -l (8) 
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We are interested in odd solution of (jHJ) which goes to =pl on infinities, i.e. 

lim <p(t) = =Fl (9) 

t— »±oo 

and 

^(t) = ^(-t) (10) 

2.1 Construction of Solution Using Iterative Proce- 
dure 

To construct solution of integral equation © with the properties ([9*|)-([T0"J) 
one could use the following iterative procedure [21 E] 

<p n+1 = (/C^,) 1 / 3 , (11) 

where zero approximation is taken as 

<p (t) = -e(t), (12) 

where e(t) is a step function defined by 

( - 1, for t < 
e(t) = < 0, when t = (13) 
{ 1, for t > 

Note that in the rhs of (|12[) expression of the forme a 1 / 3 denotes an arithmetic 
cubic root of a which is well defined for negative arguments as 



a 



Results of the iterative procedure ()li p -()12 p are presented on figHJ 

In the next section we will prove convergence of iterative procedure ([11)1- 

(jUj) that supports the results of numerical computations in 

The first approximation ipi (t) could be computed analytically and is given 

by the arithmetical cubic root of the error function. Indeed, 

^(t) = (/C^ )(t) = -erf(t), (15) 

where 



2 r 

erf(x) = —= / exp(— t )dt 
Jo 



1 



-0.5 



Figure 1: Results of iterative procedure (|TT |) - (fT^j) for a large number of steps 
(~ 10 5 ). 



In order to prove the convergence of iterative procedure (JTTjl it is conve- 
nient to take into account the parity property (|10|). This allows us to rewrite 
equation (JHJ) on the semi-axis 

(/C_<^)(t) = <^(t) 3 , t<0 (16) 



where integral operator /C_ is defined by 



K^{t,t')ip(t')dt' 



(17) 



e -(t-t') 2 _ e -{t+t') 2 



(18) 



where the kernel K_(t,t') is given by 

V 71 " 

Let us note that there is the following positiveness property of the kernel 

K_(t,t') > for all t < 0, t' < (19) 

and 

K_(t,0) = K40,t') = 0, 

see figl21 



2.2 Convergence of Iterative Procedure 

Here we describe the construction of solution for equation (|16|) on the semi- 
axis (t < 0). Consider the following iterative procedure 



ip n+ i = (/C_y9 n ) 1/3 



(20) 
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Figure 2: The kernel K_(b,t) function of t for b 
positive for t < 0, b < 0. 

Let us note, that the second iteration 

ip 2 = (KL<^i) 1/3 
could be represented in the form 

<p 2 (t) = < Pl (t)(l-A(t)), 

where A(t) is defined by 

<Pl(t) -<P2(t) 



-0.5 



-5,-0.5,-0.1. It is 



A(t) 



<Pi(t) 



(21) 



(22) 



(23) 



On figJ3]the difference ((pi(t) — <f2{t)) an d A(i) are shown for negative t. 
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Figure 3: Plots of tp^t) - (f>z(t) and A(t) illustrate that A(t) < 0.05. 
From (fT5jl . (|2T1). and pty it follows that (see. fig|3J) 

A(t) < A max = 0.05 (24) 

thus 



(25) 



Since we have the positiveness property (|T§j) we can integrate the inequality 
(f231) with the kernel K_(t,t') 



/U pU 
K^(y, x)(fi(x)dx > / K_(y, x)(p2(x)dx > (26) 
-oo J — oo 

>(1-A max ) K_(y,x)v9i(x)<ix 

./ — oo 

Inequality (pjj) leads us to 

<pUv) > <P&V) > - A m a.) (27) 

Now taking the arithmetical cubic root we get 

My) > My) > Mv)0- ~ &m ax ) 1/3 (28) 

and more over 

My) > My) > My) > MvK 1 - A™-) 173 > MvX 1 - A ma *) 1+1/3 (29) 

Analogously we have 

Mv) > My) > Mv)£ - A ma .) 1+1/3+(1/3)2+ -+ (1/3) "- 2 = (so) 

= (px(y)(l - A maa .)i~2(!) n ~ 2 ; 



i.e. 

Mv) > My) > MyX 1 - A max )iM\) n -\ (3i) 

From (|31|) it follows that f n {y) is uniformly bounded on the whole negative 
semi-axis. 

More over, beside (|27jl and (|28|l we have 

(i - A max ) 1/3 My) < vKv) < vKv) (32) 

thus 

(i - A max ) 1/9 My) < My) < My) (33) 

Analogously we have 

(i - A max ,) 1/3n_ Vn(2/) < VnMy) < My) (34) 



Hence 

\<Pn(v) ~ <Pn+l(y)\ < <fn(y) (l - (1 - A max )^^j (35) 

Finally, taking into account that (p n (y) is uniformly bounded we obtain 

C 

Wn{y) - <fn+i(y)\ < ^, (36) 

where C is constant. 

From it follows uniform convergence of (f n {y) on the semi-axis, that 
is uniform convergence of the iterative procedure on the semi-axis. 



3 Scalar Field Dynamics in Fermionic String 
Model I 

3.1 Integro-Differential Equation 

In this section we consider the following integro-differential equation 

(-g 2 ^ + l)(/C^)(t)=^(t) 3 , (37) 

where the integral operator JC is defined by (jHJ) and q is a parameter. This 
equation describes dynamics of the homogenous scalar field with the lowest 
mass square (tachyon) in fermionic string model in some approximation, 
please see section E] for more physical details. The equation (pTfjl transforms 
to p-adic equation © in the case q — 0. Although the value of parameter q 
for fermionic string model is given by 

String = "IT^ « ' 96 > ( 38 ) 
4111 3v/3 

we consider equation (f3Tj) for various values of parameter q. 

3.2 Fully Integral Form and Iterative Procedure 

The equation (|3*7j) could be written in fully integral form 

(K 9 <p){t) = <p(t) 3 , (39) 
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where the kernel of integral operator K, q is obtained by differentiating the 
kernel of (j5J), namely by applying operator (—q 2 d 2 + 1) to exp[— (t — t') 2 } 

i r°° 

{K q <p){t) = -= / e~^ 2 [l + 2q 2 (l - 2(t - t'f)] V {t')dt' (40) 

V 71 " J-oo 

As in the previous section we are interested in odd solution of (|39j) which 
goes to =Fl on infinities, i.e. the solution of (p?9|) with the properties (J9|)-([TU]). 
To construct the solution we use the following iterative procedure 

<p n+1 = (/C^) 1 / 3 , (41) 

where zero approximation is taken as 

<p (t) = -e(t), (42) 

and e(t) is a step function defined by (|T3j). In the previous section we proved 
that this iterative procedure converges to a solution for the case q = 0. 

3.3 Numerical Results 

We use the iterative procedure (|41| ) -(|42 p to perform numerical investigation 
of the solution properties for various values of parameter q. 

It was found that there exists a critical value of the parameter q 2 : q 2 r 
1.37. This is the maximum value of q for which there exist interpolating so- 
lutions. For q 2 < q 2 r there were numerically found solutions with asymptotic 
behavior Q which oscillate along —e(t) with exponentially decreasing am- 
plitude. On figHjit is demonstrated a numerical limit of iterative procedure 
@D-g2J for q = q string . 

While increasing the value of q above the critical value q cr there appears a 
swing and the solution of (|3~9|l becomes periodic (see figEJ). Let us note that 
this swing for q 2 ~ 1.37 ± 0.01 appears on quite big step numbers (values of 
n in (|4ip). On figlHlthe result of iterative procedure for q 2 = 1.38 for various 
step numbers are shown. 

Numerical investigation of iterative procedure on very high step numbers 
(~ 10 5 ) showed that transformation to periodic regime does not appear for 
q < q cr . Although for values of q a little above the critical value q cr the 
swing, i.e. a transformation to periodic regime appears on much smaller step 
numbers (~ 10 2 ). This shows that there is a principle difference between 



10 



0.5 



-4 



-2 



2 



4 



-0 . 5 




Figure 4: Results of iterative procedure (fH |) -P% J) for q = q s t r in g - 

solutions for q below and above q cr . This fact is discussed in more details in 
the next section. 

Here we used (|42|) as a zero approximation (p which is not continuous 
at the point t = 0. To understand that this discontinuity does not affect 
the results of iterative procedure we tried several smooth continuous zero 
approximations which have asymptotic behavior Q. In particular we tried 
— - arctan(t) and — erf(i). The results of the iterative procedure with these 
zero approximations where for large step numbers the same as with (|42p. 

For numerical computation we used the following approximation for K, 



Here A was automatically adjusted, namely A ~ 10. This approximation 
lead us to the algorithm with linear complexity that allowed us to compute 
~ 10 5 iterations on parallel 64-bit parallel Sun Ultra-SPARC machine. The 
entire algorithm was written in C++. 

To construct numerical algorithm we essentially used object-oriented de- 
sign [01 Elj. This allowed us to develop a general algorithm which solves 
some class of integral equations, in particular it was used to solve a system 
of nonlinear integral equations discussed in section HJ 



t+A 




(43) 
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Figure 5: Appearance of the swing of the solution with the increase of q. The 
case q 2 = 1 (thick line) - no swing, interpolating regime, q 2 = 1.4 (medium- 
thickness line) - appearance of the swing, oscillatory regime, q 2 = 1.8 (thin 
line) - oscillatory regime. 



j 




A 


I 




A 




A 


l 


h 

1 




A 




A 




A 


A 


A 


A 




















\i 




5 


























20 






L0 






-0 


5 




^ 






1 




I 


20 








li 


I/ l 


j 


V 




V 




V 


-1 
.5 


i 


1/ 




1/ 






i I 


i 


1 



Figure 6: Transformation to the periodic regime for q 2 = 1.38 appears on 
quite big step numbers - 100-th step (thick line), 250-th step (medium thick- 
ness line), 500-th step (thin line). 
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3.4 Two Regimes of The Solution 



In this section we are interested in the mechanism which forms numerically 
found exponentially decreasing oscillations and the presence of q cr . The basic 
idea is to present a solution with nonzero q as a deviation along the solution 
with q = 0. We write linear equation for this deviation for large values of t. 
Let us write the solution of (|37jl as a sum 

<p(t) = Mt)+x(t), (44) 

where (fio(t) denotes the solution of (J1J. Substituting to (|3*7jl and leaving 
only linear terms in x we get 

(-q 2 d 2 + 1)/C(0 O + X) = 4>l + Hlx (45) 

Now using the fact that 0o satisfies (j3J) we get the linear integro-differential 
equation on x(t) 

(~q 2 d 2 + l)K X = m X + q 2 d?Wo (46) 

From ()43|) we see that in K, integration could be taken in finite limits with 
good precision. This gives us ability to write large t approximation of (J4l)j) 

(-q 2 d 2 + \)K X = 3x, (47) 

here we used the fact that 4>o(t) ~ 1 and its derivatives are zero for large t. 
Representing x{t) as 

X (t) = e im 

and substituting it to (|T7|) we get the following characteristic equation 

{q 2 Q 2 + l)e^ n2 = 3 (48) 

We consider (J48|) as an equation for complex variable Q with parameter q 2 . 
There is a minimum value of q 2 : q 2 ~ 1.77 for which Q is real. For q 2 < q^ 
it has solutions with nonzero imaginary parts which gives oscillatory regime 
with exponentially decreasing amplitude. For q 2 > q 2 , solutions are real. As 
it was mentioned in the previous section numerical computations give the 
critical value q 2 r ~ 1.37 which is smaller than found above value q 2 , ~ 1.77. 
This probably means that one can not neglect the difference of 0o from 1 
in the rhs of (jH^) or even one has to take into account nonlinear terms in 
X in (J43j) . Nevertheless method discussed in this section provides a good 
qualitative explanation to the behavior of solutions for small values of q. 
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3.5 Periodic Solutions For Large q 

In this section we investigate the asymptotic behavior of the solution for large 
q. First, let us introduce a new function x(t) defined by 

X(t) = f(qt), (49) 
In terms of x{t) the equation (J37)l rewrites as 

(-d? + l)(V qX )(t) = x(t)\ (50) 
where integral operator V q is defined by 

(V qX )(t) = 4= / e-« 2 Wx{t>)dt (51) 

V 71 " J-oo 

Noting that solution of (|5()j) depends on parameter q, let us write it as a 
power series in 1/q 2 . Introducing a change of variables t' — > r = t + t'/g in 
the rhs of (joTj) we get 

q 



/GO 
e-" 2r2 x(t + r)rfr 
-co 



'7T 

Now changing variables again r — > a = gr we get 

1 /"CO _ 

1 / „2 . a 



e a x(t + -)rfcr 

V-oc 9 



Finally, expanding x(t + cr/q) in Taylor series with respect to o we get 

-F= / e-°\ X (t) + -X'(t) + -X"(t) + ■ ■ -)da 
J-oo Q T 

All terms which are the odd powers of a vanish after integration. We obtain 



or 

d 2 



__ + l) e 3?^ x = x 3 ; (52) 
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where exp(l/4g 2 <9 2 ) is understood as a formal expansion 

3^ 9 * = V ( — — 

Let us note that (|52|) could be directly obtained from differential form of the 
equation, see sectional 

In a limit q — > oo the equation (|52jl becomes an equation for anharmonic 
oscillator. This explains numerically found periodic solutions for large q. 



4 Scalar Field Dynamics in Fermionic String 
Model II 

4.1 System of Nonlinear Integral Equations 

In this section we study the following system of integral equations 

(K*){t) = if(tf (53) 
(-g 2 ,9 2 + l)(/C^)(t) = a(t)<p{t) (54) 

where the integral operator K, is defined by (JSJ) and q is a parameter. Phys- 
ically we are interested in solutions which are finite on the whole axis and 
<p(t) is odd function interpolating between =)=1. This system of equations de- 
scribes dynamics of the homogenous scalar field with the lowest mass square 
(tachyon) in fermionic string model in the first nontrivial approximation, 
please see section El for more physical details. 

We are interested in the solutions with the following properties. The 
function (p(t) is odd and goes to =)=1 on infinities, i.e. satisfies (JUJ-tfTUJ) and 
the function a(t) is even and goes to 1 on infinities, i.e. 

a{t) = a(-t), a(Too) = 1 (55) 



4.2 Reduction to a Single Equation and Iterative Pro- 
cedure 

In order to build the solution of the system (|53| ) -([54 p first we rewrite it as a 
single nonlinear integral equation 

;-g 2 d 2 + w 



<P Z = 1C 



(56) 



15 



This allows us to construct the following iterative procedure 



¥>n+i(*) = ~£{t)> 



(57) 



where integral operator fC q is defined by PO]). and zero approximation is 
taken as 

<p (t) = ~<t), (58) 

where e(t) is a step function defined by (|T3|). Please note that although here 
we discuss iterative procedure for <^(i) one can also write a similar procedure 
for a(t) taking as a zero approximation a continuous function o"o(t) = 1. One 
could see that (|57|) assumes that </? n ^ that makes us define y?o at t = in 
a special way, see next section for further details. 



4.3 Numerical Results 

We used iterative procedure (joT|) to investigate solutions of the system 

As it was mentioned in the previous section we used — e(t) as a zero 
iteration. Since in (J57|) we have to inverse f n to build the (n + 1) iteration 
we define in the rhs of ()57|) and (J58|) to be nonzero for t = 0. The 
basic idea here is to make the value of e(0) to randomly 1 take values ±1. In 
particular we used the following relation 

e(0) = (-If (59) 

where n denotes the iteration number. We also tested some more complex 
random sequences, but iterative procedure with (jSUjl was the fastest and 
more predictable. One could also take e at t = to be constant for example 
e(0) = 1 but this lead to some asymmetry in the resulting solutions. 

The results of iterative procedure (|57ji for q = q s trin g are presented on 
figd Please note that the solution (p has a break at t — 0. 

It was found that there exists a critical value of the parameter q 2 : q 2 cr pa 
2.24. This is the maximum value of q for which there exist interpolating 
solutions. For q 2 > q 2 r the iteration procedure starting from some step faces 
the negative argument under the square root in the rhs of (|57j). This means 

1 The author is grateful to L.V. Joukovskaya for the idea to use randomness in deter- 
ministic algorithms. 
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Figure 7: Results of iterative procedure solving the system fit) 
(thin line) and a(t) (thick line) for q 2 = 0.96. 

that iterative procedure (|H7j) while valid for finding interpolating solutions 
fails to find oscillatory ones. 



4.4 Linearization of the System for Large t 

In this section we investigate the behavior of the system (|53 )) - (J54|) for q < q cr 
in the large t limit. The basic idea is analogous to section EOl We represent 
the function tp which forms the solution for q ^ as a deviation from (p - 
the solution for q — 0. We write a linear integral equation for this deviation 
in the large t limit. 

Let us write the solution of (JHHjl - lJHl)) as a sum 

¥>(*) = ¥>o(t)+x(*), (60) 

where ifo(t) is the solution for q = 0. Substituting to (|55)i and leaving 
only linear terms in x(t) we get 

AC— 2 — + /C^ AC X^2~ = <^o + 2l PoX (61) 

Using the fact that <^ solves the system for q = and <^o(*) — 1 and its 
derivatives are zero for large t we obtain 

/CAC gX - ACx = 2x (62) 

Representing x(t) as 

X (t) = e** 
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and substituting it to (|62|) we get the following characteristic equation 



(q 2 n 2 + l)e 




2 



(63) 



Analogously to what we did in section 13.41 we consider (JfiHj) as an equation 
for a complex variable Q with parameter q 2 . There is a minimum value of 
q 2 : q 2 ) « 3.05 for which Q is real. For q 2 < q 2 it has solutions with nonzero 
imaginary parts which gives oscillatory regime with exponentially decreasing 
amplitude. For q 2 > q 2 « 3.05 solutions are real. As it was mentioned in the 
previous section numerical computations give the critical value q 2 r m 2.24. 
This difference probably means that one can not put ipo equal to 1 in the rhs 
of (Jfi2|) or even consider nonlinear terms in \ i n ((SI]) - 

4.5 Asymptotic Behavior For Large q 

Here we investigate asymptotic behavior of the system (|53 )) - ()54)) for large q. 
As in the section 13.51 we introduce a change of variables analogous to (|4*H|) 



where integral operator V q is defined by (|5ip. 

Performing computations analogous to section !3.5l we get that in the large 
q limit the system (jMJ)-(jnSJ) becomes 



x(t) = f(qt) 




x{tf 



(64) 
(65) 



{-dl + l) X (t) 



x(tf 
ttt)x(t) 



(66) 
(67) 



that is equivalent to an anharmonic oscillator. 
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5 Differential Form of the Equations and Phys- 
ical Roots 

In this section we provide physical details including actions for which corre- 
sponding equations of motion form the subject of this paper. 
Effective p-adic action is given by (21 El H] 



S 4 /' <'"■<■ 

9p 



-(pp ^(j) 

2 p 



1 p 2 

g'p 9 2 p-^ 



,- , , (68) 



Here is a scalar field which describes tachyon in p-adic string model, x = 
(t, af ) are d-dimensional space-time coordinates, p is a prime number, g p is a 
coupling constant (g is universal coupling constant), the D'Alamber operator 
is defined in a standard way 

□ = -— + V-V, (69) 
and operator p~^ D is understood in the sense of expansion 

n=0 ^ ' 

The corresponding equations of motion are the following 

p-l n ( f) = ( f)P (71) 
For homogenous field configurations (|7Tj) rewrites as follows 

v ¥1$ = <f (72) 

If we consider slowly varying solutions and neglect high order derivatives 
in the lhs of (J72j) we get equations for an anharmonic oscillator 

l -\Yipd 2 t <P + ci> = <f, (73) 



in the potential 



™ = £ & - ?TT^ +I ' (74) 
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Figure 8: p-adic potential for the case p = 2 (left) and p = 3 (right). 



From (J74"j) we see that the cases with even and odd p we have qualitatively 
different behaviors jH]- 

In particular in the case p = 2 the potential (|74|) has minimum when 
= (fig HI) and maximum when = 1. In jHj it was proved that there are 
no monotonic solutions interpolating between these points. 

In the case p = 3 the potential has three extremal points 



(,(1,3) 



o. 



4 3) 



1 



(75) 



The extremal points 0o 1 ' o; correspond to unstable vacua, and 0q 2 ^ corresponds 
to the stable one. In the previous sections we were interested in the time 
dependant solutions symmetrically interpolating between 0^ and 0q 3 ' ) . 

The action for the scalar field with lowest mass square (tachyon) in 
fermionic string model in the first nontrivial approximation is given by [S] 



S[v,1>] 



dt 



where 



and 



T(t) = exp(--d 2 )v(t), (t) = exp(--9 2 )^(t) 



a 2 

^string 



4 In 7 



0.96 



(76) 



(77) 



Here ip is the tachyon filed and v is the auxiliary filed. We will write the 
equations of motion in terms of T and \I/ while the physical fields ip and v 
will be obtained using the inverse of (J77J). 

First let us consider an approximation for the action (fTtjj) which leads to 
simpler equations of motion. If we assume that we can neglect the smoothness 
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of auxiliary field, i.e. in the interacting term write v instead of T we get the 
following approximate action [S] 



S[v,ip] = / dt 



l v{tf + ^(d t m) 2 + U 2 (t) - lv(t)* 2 (t) 



(78) 



4 w 2 v T w ' 2 r w 2 
This action leads to the following equations of motion 

v(t) = * 2 (t) (79) 
{-q 2 d 2 + l)e* 9 ?V{t) = v(t)V(t) (80) 

The first equation gives the expression for v in terms of substituting it to 
the second one we get 

(-q 2 d 2 + l)e^tf(t) = ^(t) 3 (81) 

If we rewrite f!81|) in the integral form (see (J3J) with a = 1/4) we get 

(-q 2 d 2 + 1)/Ctf = ^ 3 (82) 

This equation was studied in section El see equation (J37j) . 

The equations of motion for original action (|76J) are the following 

e^ 2 T(t) = ^ 2 (t) (83) 
(-q 2 8 2 + l)e^(t) = T(t)*(t), (84) 

In the integral form this system is equivalent to the system ([53 )1 - ([54 )1 which 
was studied in section EJ 

Let us compare the dynamics of physical field ip(t) obtained from exact 
action (|75jl and approximate action (|75|) . We compute physical field 
from solutions of ()81|) and ()83 |1 -()84 |1 using the inverse of (J77J) which has the 
following form 



^(t) = e^^it) = y°° e-^-'') 2 ^')^' (85) 

Note that although ^ solving (J8*3jl - (j8*3)l has a break at t — (see section IPj) 
the application of (|53jl makes the resulting physical filed ip smooth, see figHH 
The results of comparison of physical fields obtained from exact and ap- 
proximate actions are presented on figlTUl 



21 



Figure 9: Physical fields ip(t) (thin line) and v(t) (thick line) are smooth for 

all t (q 2 = qltring)- 




Figure 10: Comparison of physical field dynamics ip(t) obtained from exact 
action (fTHJ) . thick line, and approximate action (J7%|) . thin line, for q 2 = q 2 tring - 

Conclusions 

In this paper we numerically studied nonlinear equations depending on one 
variable with infinite number of derivatives. Two different regimes for solu- 
tion of equations of motion for the fermionic string model are found and the 
critical value of the parameter is derived. It would be interesting to explore 
in more detail the transition between two regimes. Also generalization of the 
results of the paper to the case of several variables would be important. 
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